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On 

^ ! Abstract 

Measurements of dilepton production from charm decay and Drell-Yan processes 
Q^ ■ respectively probe the gluon and sea quark distributions in hadronic collisions. In 

^ ■ nucleus-nucleus collisions, these hard scattering processes constitute a 'background' 

to thermal contributions from the hot matter produced by the collision. To de- 
Ph. termine the magnitude and behavior of this background, we calculate the hard 

scattering contribution to dilepton production in nuclear collisions at RHIC and 
LHC at next to leading order in perturbative QCD. Invariant mass, rapidity and 
transverse momentum distributions are presented. We compare these results to op- 
timistic hydrodynamic estimates of the thermal dilepton production. We find that 
charm production from hard scattering is by far the dominant contribution. Ex- 
periments therefore can measure the gluon distribution in the nuclear target and 
projectile and, consequently, can provide new information on gluon shadowing. We 
then illustrate how experimental cuts on the rapidity gap between the leptons can 
aid in reducing the charm background, thereby enhancing thermal information. 
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1 Introduction 



Dilepton production provides an important tool for measuring the temperature of the 
high density matter produced in the early stages of a relativistic heavy-ion collision (see 
|l[ and references therein). To make use of this tool at ion colliders such as RHIC and 
LHC, we must understand additional sources of dileptons from hard scattering and other 
non-equilibrium processes. In high energy pp interactions, the continuum in the dilepton 
mass range above M ~ 2 GeV is dominated by the Drell-Yan process and by semileptonic 
decays of D and other charm mesons. These leptons are produced by hard scatterings, 
at scales exceeding M and 2m c ~ 3 GeV respectively, so that their production can be 
addressed using perturbative QCD. In contrast, the lower mass region arises from soft 
processes, for which theory offers little guidance. We therefore hope to find a signal of 
thermal dilepton production at masses M ~ 2 — 3 GeV where i) the thermal contribution 
can still be sizable and ii) the background is calculable. 

In this paper, we compare the predictions of dilepton production from a simple ther- 
mal model with the hard 'background' from heavy quark pair, QQ, decays and Drell-Yan 
production at RHIC and LHC nucleus-nucleus collision energies, yfs = 200 GeV and 5.5 
TeV in the nucleon-nucleon center of mass. Previously we compared the thermal dilepton 
and thermal charm rapidity distributions with Drell-Yan and initial charm production [p]] 
calculated at leading order (LO). Our new next-to- leading order (NLO) results, shown to 
agree with pp and pp data in H|, no longer depend on arbitrary phenomenological 'K 
factors' (often incorrectly taken to be '2' in the literature). The remaining uncertainties 
in the NLO perturbative approach are well defined and likely rather small 0. Also new 
to this paper are calculations of the invariant mass and transverse momentum distribu- 
tions of the lepton pairs. We use these distributions to illustrate how thermal and hard 
dileptons can be distinguished in an experiment. In addition, we introduce the dilepton 
contribution from B meson decays, and discuss the effects of nuclear shadowing on the 
initial production. 

Figure 1 shows the invariant mass distributions of the calculated contributions to the 
dilepton continuum in central nucleus-nucleus collisions at RHIC and LHC. We find that 
hard charm quark production and decay dominates the continuum below the T mass. 
In particular, the charm signal is more than an order of magnitude above the optimistic 
thermal dilepton and thermal charm rates for M > 2 GeV. This result implies that 
dilepton measurements can be used to extract the low x gluon density in the nucleus (see 
also ||). On the other hand, the isolation of thermal signals will not be straightforward. 

Additional complications can arise from the fact that the charm production cross 
section is large enough for multiple DD pairs to be produced in a single nucleus-nucleus 
collision. Uncorrelated pairs form when a lepton l + from one DD is randomly paired 
with a l~ from another DD; the correlated signal in Fig. 1 includes only dileptons from 
DD pairs in which both quarks decay to leptons. In central Au+Au collisions at RHIC, 
up to 67 uncorrelated pairs may be produced while at the LHC, over 3000 uncorrelated 
pairs will contribute to the raw continuum from Pb+Pb collisions. Uncorrelated charm 
does not affect the measurement of parton densities, but it further complicates the task 
of extracting thermal information. Ideally the uncorrelated pairs can be removed by a 
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like-sign subtraction, leaving only the pairs shown in Fig. 1. Whether this subtraction 
can hold at the accuracy needed to extract thermal signals is another matter. 

To enhance the thermal signal in an experiment, one can choose to count only those 
lepton pairs that have a small separation in rapidity. As observed by Fischer and Geist, 
dilepton pairs from charm decays typically occur with a large rapidity gap [f| . No such gap 
is present in Drell-Yan pairs or - importantly - in thermal pairs. Together with like sign 
subtraction, we expect that rapidity gap cuts can essentially remove the uncorrelated 
charm contamination and greatly suppress the correlated charm background. Indeed, 
the finite acceptance of a real detector can serve a similar purpose. We show that the 
acceptance window of PHENIX/RHIC and ALICE/LHC can enhance the signal from 
thermal charm decays to the point of measurability by rejecting pairs with large gaps. 

This paper is organized as follows: In section 2, we discuss in detail Drell-Yan, DD, 
and BB production in pp interactions and comment on how these results may be modified 
by nuclear shadowing. In section 3, we show the results for thermal dilepton and thermal 
DD pair production at RHIC and LHC, assuming the most optimistic scenario for the 
initial conditions to maximize the thermal rate. We compare the initial hard production 
of dileptons to the thermal model results in section 4. We discuss the effects of realistic 
detector geometries on the dilepton spectra at RHIC and LHC and draw our conclusions. 



2 Dilepton Production from Initial Interactions 

Perturbative QCD calculations of Drell-Yan and heavy quark production at leading or- 
der have long been available. The LO calculations differed from the experimental mea- 
surements by a i^-factor (K = a cxp / a thcory ) of 2-3 for charm production and 1.5-2 for 
Drell-Yan production. This difference suggested that higher order corrections to the pro- 
duction cross sections were important. Additionally, while multiplying the leading order 
cross section by a K factor describes the single-inclusive quark distributions as well as 
the mass and rapidity distributions of Drell-Yan and QQ pairs, the pair px distributions 
for both processes are trivial at leading order since the pairs are produced back-to-back. 
Therefore next-to-leading order calculations are necessary to fully describe hard dilepton 
production. 

In our previous work we used the LO cross sections. Since then, an NLO treatment 
of QQ production has been made available ||. An NLO treatment of the Drell-Yan p? 
distribution is also now available || [7|]. With the NLO description of the perturbative cross 
sections, we also use NLO evaluations of the parton densities. Our results are obtained 
using the MRS D— ' || parton densities^, compatible with the low x data from HERA 
[Tofl . The MRS D— ' sea quark and gluon distributions grow oc x~ 1//2 at the initial scale, 



Qq = 5 GeV 2 , when x — > 0. Using recent parton distribution functions that agree with 
the HERA data produces a substantial increase over our previous results JT|, obtained 
with obsolete leading order parton distributions that become constant at Ql as x — > 0. 



Observe that older parton distributions such as [i t] substantially underestimate the initial 



1 All available parton distribution functions are contained in the package PDFLIB ||, available in the 
CERN library routines. 
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Drell-Yan and QQ production at heavy-ion colliders. 

We remark that updated versions of the MRS distributions are available, including 
MRS G [0 which has a slower low x growth than MRS D— '. While small, the changes 
from D— ' to G affect our charm results most strongly. Compared to D— ', we find that 
charm rates from G are 5% larger at RHIC and 20% smaller at LHC. No doubt, these 
numbers will continue to improve at that level as more data is analyzed, both from HERA 
and from Fermilab |p"3f . 

While the NLO evolution generally improves the agreement between the theory and the 
data, additional uncertainties are introduced, including dependence on scale and scheme. 
At fixed order in perturbation theory, the calculations depend on the renormalization 
scale and the factorization scale as well as the regularization scheme (MS or DIS). The 
hard scattering matrix elements and the definition of the parton densities for each process 
are specified by the regularization scheme. In this paper, we have used the MS scheme. 
The renormalization scale enters into the strong coupling constant, a s , and the partonic 
cross sections while the parton distributions are evaluated at the factorization scale. The 
precise relationship between these scales and the momentum transfer, Q, is not uniquely 
defined. However, since the parton densities are analyzed assuming that they are equiv- 
alent, we also assume this. If the perturbative expansion converges, further higher-order 
corrections are small at large enough values of the scale \i. For such scales the phys- 
ical cross section should become independent of the scale and scheme when calculated 
at higher and higher orders. If the /i dependence is strong, the perturbative calculation 
at that order is unreliable and further higher-order corrections are necessary [FJ]. The 
rates from the initial hard scatterings are rather sensitive to the scale and scheme in the 
moderate mass and regime relevant to heavy-ion experiments. 

Another uncertainty in our results involves the nuclear dependence of hard processes. 
When the charged parton distributions are probed in deep-inelastic scattering with a 
nuclear target and compared to a deuterium target, the ratio Rf 2 = F^/F® has a char- 
acteristic shape as a function of x. A depletion in the nucleus is observed at low x, the 
shadowing region, and intermediate x, the EMC region. Shadowing occurs in the region 
below x ~ 0.1 while the EMC region refers to 0.3 < x < 0.7. Between the shadowing 
and EMC regions, Rf 2 > 1, referred to as antishadowing. Although the origin of this 
behavior is not well understood, it is postulated to be either an interplay of coherent 
and incoherent multiple scatterings in the target or a modification of the parton densi- 
ties in nuclear matter. In any case, the effect can be modeled phenomenologically by a 
parameterization to fit the nuclear deep-inelastic scattering data and implemented by a 
modification of the parton distributions in the nucleus. In the appendix, we choose two 
different parameterizations of the nuclear parton densities to illustrate the effect in AA 
collisions at RHIC and LHC. The first is a general fit to the most recent nuclear target 
data |TJ|] that does not differentiate between quark, antiquark, and gluon modifications 
and does not include evolution in Q 2 . The second modifies the valence and sea quark 
and gluon distributions separately and includes Q 2 evolution |TB[ but is based on a fit to 
somewhat older data using obsolete parton densities. We find that at RHIC energies the 
charm and Drell-Yan yields may be reduced 40-50% while the bottom yield is changed 
by « 10%. At the LHC, all the yields are reduced 50-60%. The results depend on the x 
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region probed and, in the case of the second parameterization, the scale Q 2 . 



2.1 Drell-Yan Production 



A detailed discussion of Drell-Yan production in high-energy NN collisions can be found 
in Ref. ||. We repeat some of the pertinent points here. At leading order, Drell-Yan 
pairs are produced by qq annihilation into a virtual photon which decays to a lepton pair, 
QQ ~~ * 7* ~~ * l + l~- The production cross section for lepton pairs with invariant mass, M, 
and rapidity y, summed over quark flavor / is 



M 



da 



DY 



dydM 2 



^or^e)[q f (x l ^)q f (x 2 ^) + (1 -> 2)] 



(1) 



where Bq = 47ra 2 /9M 2 is the LO qq annihilation cross section, r = M 2 /s, \i is the scale, 
and xx t 2 = \pT£ ±v are the projectile and target momentum fractions at which the parton 
densities are evaluated. For y > 0, xi increases with rapidity while x 2 decreases. At NLO, 
the Compton and annihilation processes qg q'j* and qq — > (77* contribute in addition 
to vertex corrections to the LO cross section so that eq. (1) is replaced with 
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^2e 2 f [q f {x u fi)q f (x 2 , fi) + (1 -> 2)] 
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S(l-z) + ^f q (z) 



J2e 2 f [g(x!, n)(q f (x 2 , n) + q f (x 2 , fi)) + (1 -> 2)] 
. / 



2tt 



Note that going to NLO requires a redefinition of x± and x 2 since the unobserved parton 
contributes to the total momentum in the final state. The correction terms f q and f g 
are regularization scheme dependent. The scale and scheme dependences are not large 
for M > 4 GeV. We will only consider Drell-Yan production for M > 2 GeV since 
below this value the perturbative calculation becomes unreliable. The mass and rapidity 
distributions were calculated using a program provided by Rijken and van Neerven [17] . 

The pt dependence is trivial at LO — the lepton pair has pt = if no intrinsic parton 
Pt is included. Some of the Drell-Yan pair transverse momentum can be accounted for 
by introducing a soft intrinsic px distribution with (p?) ~ 0.3 GeV. However, this is too 
small to account for the measured px distributions with (p T ) ~ 1 GeV at M = 10 GeV. 
At NLO, the pair acquires transverse momentum through recoil of the virtual photon 
with the final-state parton. Here the pt dependence can be calculated perturbatively and 
is well behaved for pt ~ M. However, at low pt, the perturbative expansion parameter, 
a s ln 2 (M 2 /p^), becomes large and the expansion breaks down, making it necessary to 
resum the perturbation series. The resulting cross section is 



M 2 



da r , 



dp 2 T dydM 2 



na Te, 



d 2 b 



W(b) 



(3) 
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where a Fourier transformation is made to impact parameter space and the form factor, 
W(b), effectively sums the leading and sub-leading logarithms. The pt distributions have 
been calculated using a code developed by Arnold and Kauffman for pp — > W,Z ]F] and 
extended to Drell Yan in pp collisions in [fjj . This code includes a method of interpolating 
between the low and high px behavior. 

The number of Drell- Yan pairs of mass M produced in central nucleus-nucleus colli- 
sions, 

dN DY dfjpY 

W" Tab(0) W (4) 

is always small in the range where the perturbative calculation is expected to be applicable, 
i.e. M > 2 GeV. However, only a calculation of the mass-integrated cross section can 
reveal whether or not the AB rate is large enough for false Drell- Yan pair production to be 
a problem. Since such a calculation is not possible over all masses, experimental techniques 
must be used to insure that these uncorrelated Drell- Yan pairs do not contaminate the 
spectrum. For example, a pt cutoff should reduce the probability of false high mass 
Drell- Yan pairs arising from low mass pairs with large rapidity separation. Additionally, 
a study of the angular distribution would also reduce any uncorrelated Drell- Yan pairs. 
The number of lepton pairs in central Au+Au collisions at M = 2, 4, and 6 GeV are given 
in Table 1 for RHIC and in Table 2 for central Pb+Pb collisions at the LHC. They are 
significantly larger than in our previous work |l] due to the different small x behavior of 
the parton densities. The average mass and pt of the Drell- Yan pairs (with M > 2 GeV) 
produced at RHIC and LHC energies are shown in Table 3. 

Figures 2 and 3 show the NLO rates for Au+Au collisions at RHIC and Pb+Pb 
collisions at the LHC. The calculated pt distributions for M = 4 and 6 GeV at y = 
and 2 are shown in Figs. 2(a) and 3(a). The distributions at higher rapidity are nearly 
parallel until the kinematic limit is approached. Given this, we have assumed that the pt 
distribution has the same shape for M = 2 GeV and scaled the 4 GeV results to obtain 
the pt distribution at M = 2 GeV. The rapidity distributions are shown in Fig. 2(c) 
for RHIC energies and 3(c) for the LHC. The cross section grows with rapidity until the 
kinematic limit is approached because x 2 g(x2,/u) increases and xiq(xi, fx) decreases with 
increasing y. As previously discussed, the parton density dependence is enhanced for low 
masses and high rapidities at = 5.5 TeV even though the rapidity-integrated results do 
not depend strongly on the parton densities. Note that the strongest dependence on the 
parton densities occurs outside the range of current measurements. In our previous work, 
the variation with rapidity was weaker because we used obsolete parton distributions 
where xq(x — > 0, Qq) — > constant. 



2.2 Heavy Quark Production and Decay 

At RHIC and LHC energies, heavy quark production will be substantial. We consider both 
charm and bottom production and fragmentation into heavy mesons which subsequently 
decay to lepton pairs. The double differential heavy meson pair production cross section 



6 



at LO is m 



H d 3 p H d 3 p w 



s dx x dx 2 t ^EhEjj 



27T X\ X2 

Dq/h{zh)Dq/h{,z 1 j 

Z H Z J[ 



E0E7 



$ 4 (Pi+P2-Pq-Pq) 



(5) 



where H(Qq)H(Qq) is the heavy meson pair, DD for Q = c and 55 for Q = b. The 
fragmentation function, Dq/h(z), describes the hadronization of the heavy quarks where 
z = Ph/pq is the fractional momentum of the heavy quark carried by the hadron. Charm 
hadroproduction at low px is best described by the assumption that the charmed quark 
experiences no momentum loss during hadronization, i.e. Dq/h(z) = 5(1 — z) ||18|| , result- 
ing in a shift between jjq and yn due to the quark and meson mass difference. At LO 
heavy quarks are produced by gluon fusion, gg — > QQ, and quark- ant iquark annihilation, 
qq — > QQ. To LO, the convolution function C(xi,x 2 ) is 

C(xi,x 2 ) = ^[xiqfixt, jj,)x 2 qf(x 2 , p) + x{q f (xi, p)x 2 q f (x 2 , p)} da ^ Xl '^ 2 ' | g _ (6) 

1 x / .d$(a:i,a:2,A*)| 
+ x^Xi, n)x 2 g{x 2 , p) -~ | M , 

where da/dt can be found in e.g. || [L4|]. 

The NLO, order a 3 , corrections to QQ production have been calculated for the total 
cross section f20fl , single-inclusive quark distributions ^> HU) an d exclusive QQ pair 
distributions ||. In addition to real and virtual corrections to the LO processes, quark- 
gluon scattering, q(q)g — > QQq(q) is also included^. The total partonic cross section, oV,-, 
can be expressed as 



5 : ii (xiX 2 s,mQ,/i) 



rrin 47r 



flAp) + fUp)H^/ra 



Q> 



0(al 



(7) 



where p = 4mq/ x\x 2 s and /° 
section to NLO is then 



0. The double differential meson pair production cross 



ErEjj 



da HH 



d 3 p H d 3 p 1T 



/ dxidx 2 dzndz-jj 



E H E li D Q / H (z H )D- 



Q/H\ Z H! 



EqEq 



EqEq 
da ij (xiX 2 s,mQ,p l ) 



(8) 



d 3 p H d 3 p w 



qi(xi, p)qj(x 2 , p) 



where qi and qj are the quark, antiquark and gluon densities, appropriately defined in 
the particular scheme. The sum runs over gg fusion, qq annihilation and q(q)g scattering. 
The NLO corrections become large when itiq/^/s <C 1 since gluon exchange dominates 

2 This process has been interpreted at LO as the scattering of a heavy quark excited from the nucleon 
sea with a light quark or gluon, e.g. gQ — ► gQ, and is referred to as flavor excitation (f4| po| . However, 
for moderate pt, flavor excitation is suppressed at LO due to the small heavy quark parton density near 
threshold. 
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the asymptotic behavior. However, the perturbative expansion may still be valid if further 
higher-order corrections are small. 

The cc total cross section at NLO, a^, has been compared with pp and pA data at 
*/s < 63 GeV [^3|, p4| , assuming a linear nuclear dependence |25|, |(| ^7|, to fix m c and 



fi and provide an extrapolation to collider energies ||. Reasonable agreement with the 
data was found for m c = 1.2 GeV and /i = 2m c for MRS D— ', leading to er^f = 344 /xb at 
RHIC and 17.7 mb at the LHC We use m;, = ji = 4.75 GeV to calculate bb production, 
finding = 1.5 /ib at RHIC and 224 /zb at the LHC. The cc cross sections are larger than 
earlier estimates Jl| because of the low x behavior of the gluon distributions^. However, the 
theoretical K factor, ^q^ / "^ ~ 2 — 3, is rather large, particularly for cc, suggesting that 
the perturbative expansion is unreliable for itiq/^/s <C 1. For details of the calculation 
and the theoretical uncertainties, see [|], |9) . 

It was recently shown that when /i oc m<r, the theoretical K factor is nearly constant for 
distributions that are nontrivial at LO [290. Therefore, the NLO calculation is essential 



only for the pair px and azimuthal distributions. We calculate the meson pair invariant 
mass distributions and the double differential rapidity and rapidity gap, y g = yu — Ujji 
distributions da H jj/dydM and da H jj/ dy g dM at LO assuming Dq/ h (z) — 8(1 — z) and 
checked that our results agreed with the LO results from the program of Nason and col- 
laborators ||. The rapidity gap is needed for acceptance studies since the mesons decay 
independently. The total cross section and QQ pair px distributions are calculated at 
NLO using the same program |J. We note that, so far, no resummation of the QQ px 
distribution has been performed for low px, analogous to Drell-Yan production (see eq. 
(3) and 0,0). Therefore, we have only shown the QQ px distribution where the calcu- 
lation may be considered to be reliable. The QQ and HH px distributions are assumed 
to be equal, in keeping with the trivial fragmentation function used for the longitudinal 



momentum distribution Ul8| . At large px, this is no longer a good approximation. How- 
ever, in the low px region, of interest in heavy-ion collisions, this assumption should be 
reasonable. We assume no intrinsic light parton px in the parton densities and integrate 
over X\ and 22 using four-momentum conservation. The triple differential cross section, 
eq. (5) is then 

da HlI (T NLO f 2 . > . A _n ( ,E h E 7T Dq /h (z h )Dq /1 j(z 1t ) 

' dpx^QdzHdzjfdyndyjfC (xi,X2) 



dydy g dM* a LO J 2 ' w x ' " E q Eq z hZ7T 

8{y - (y H + y7r )/2)5(M 2 - (p H + pjj) 2 )5{y g -y H + m ) , (9) 

where ctnlo/clo is the theoretical K factor. The definition of the pair rapidity, y = 
(y H + yjj)/2, holds when px,Q = p T q which is true when the initial parton p T is neglected. 
If the average number of HH pairs, 

% = Tab(0)ct H h , (10) 



3 Choosing another parton density with a lower initial scale but similar small x behavior [g8| leads to 
a somewhat smaller charm cross section at LHC, (7*2* = 6.7 mb while the cc cross section at RHIC is 351 
/ib, nearly independent of the choice of parton densities |J. 

4 Ref. J3(| reaches a different conclusion with /i oc toq. However, this choice introduces large logarithms 
for px > rriQ, making the calculation unstable. 
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produced in central collisions is small, N h jj <C 1, the lepton pairs will be correlated 
with N£° TT = N h jjB r (H/H — ► l^X). However, if N h jj > 1, opposite sign lepton pairs 
from uncorrelated HH pair decays need to be taken into account. When N H jj 3> 1, the 
average number of uncorrelated lepton pairs is N^ ncorr = N h -jj(N h jj—1)B r (HH — > / ± X). 
If Njjjj !=a 1, a distribution in N H -jj must be considered to calculate the uncorrelated pairs. 

The cc production cross sections are large enough for lepton pair production from 
uncorrelated DD decays to be substantial in nuclear collisions, particularly at the LHC. 
Given our values of cr^f and assuming that all cc pairs produce final-state DD pairs, we 
find N D p ~ 8.7 at RHIC and N D ^ ~ 450 at the LHCQ. Up to 3100 uncorrelated pairs can 
be formed among the B R (D — > l + X)N D j) ~ 55 l + and l~ leptons from the charm decays 
in central Pb+Pb collisions at the LHC. 

In Figs. 4 and 5 we show the rates of cc and DD production at RHIC and LHC 
as a function of px, M, y g , and y. The bb and BB rates are given in Figs. 6 and 7. 
The meson pair rapidity distributions are somewhat narrower than the corresponding 
quark pair distributions due to the quark and mason mass differences. The rapidity 
gap increases with M and dN/dy g dM has a minimum at y g = since for fixed mass, 
M 2 = 2m T (l + coshy 9 ), small y g corresponds to large Pt- However, the mass integrated 
rapidity gap distribution, dN/dy g , does not have a minimum at y g = 0. 

Lepton pairs from DD and BB decays have been obtained using a Monte Carlo code 
based on D decays measured by Mark-Ill and B decays observed with CLEO Op- 
posite sign lepton pairs can be produced from a single B decay by the chain B — > Dl + X, 
D — > l~X. Our calculation assumes that the lepton pairs result from the decay of the 
initial B and B. Note that the CLEO analysis assumes that the measured leptons 
are from the decays of the initial -B's rather than from secondary decays. However, if BB 
production is significant, the secondary decays in the chain could add a component to 
the low mass background. The inclusive branching ratio for D meson decay to leptons, 
averaged over charged and neutral -D's, is B R (D° / D + — ► l + X) ~ 12% |34|]. The corre- 
sponding branching ratio for B mesons of unspecified charge is Br(B — > l + X) ~ 12% 
f34|j . The momentum vectors of each meson are computed in the HH pair rest frame 
using the rapidity gap distribution, dN/dy g dM, to separate the mesons. The decays are 
calculated in the meson rest frame according to the measured lepton momentum distri- 
butions |J2[ and boosted back to the nucleon-nucleon center of mass frame where the 
lepton pair quantities are computed. To account for uncorrelated lepton pair production 
in the Monte Carlo code, two HH pairs are generated and the H from the first pair is 
decayed with the H from the second. 

The lepton pairs from correlated and uncorrelated DD decays are shown for central 
Au+Au collisions at RHIC and Pb+Pb collisions at LHC in Figs. 4 and 5. The number 
of lepton pairs from heavy quark decays at M — 2, 4, and 6 GeV are given in Table 1 

5 The cc production rate at LHC is compatible with an earlier estimate using the parton cascade model 
fnl where ~ 400 cc pairs are produced at y/s — 6.3 TeV. However, our RHIC rate is considerably smaller 
than the 60 cc pairs predicted by the parton cascade. The discrepancy cannot be fully accounted for 
within the theoretical uncertainties. In fact, if we use the same parton densities as in the parton cascade 
model with the same charmed quark mass and scale, our production rate drops by a factor of two at 
the LHC but is unchanged at RHIC energies ||. Therefore, initial cc production cannot be the only 
important source of charm in the parton cascade model. 
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for central Au+Au collisions at RHIC and in Table 2 for central Pb+Pb collisions at the 
LHC. Note that in obtaining the number of pairs from our Monte Carlo simulation, we 
average over mass bins of width 400 MeV. The average mass and pr of the decay pairs 
produced at RHIC and LHC energies are shown in Table 3. The dileptons resulting from 
uncorrelated DD pairs have larger masses since the rapidity gap between the uncorrelated 
mesons is larger on average than between the correlated pairs. In fact the average mass 
of lepton pairs from uncorrelated DD decays is 30% larger than that of the correlated 
pairs at LHC. The average rapidity of the lepton pairs from uncorrelated DD decays 
is smaller than from correlated decays at the LHC whereas the average pt is somewhat 
larger for the uncorrelated pairs. Uncorrelated BB decays have not been considered here 
but should produce only a relatively small enhancement at LHC for M < 10 GeV. We 
show the lepton pairs from correlated BB decays in Figs. 6 and 7. Note that, in general, 
the average lepton pair p? increases with dilepton mass while the average lepton pair 
rapidity decreases with mass. 



3 Dilepton Production by Thermal Processes 

We now calculate the rates of thermal lepton pair production, both directly from thermal 
qq annihilation and from the decays of DD pairs produced thermally. We expect that the 
longitudinal expansion approximately follows the scaling law, v z — z/t, at RHIC and LHC 
until large rapidities where particle densities become small and pressure gradients cause 



scaling violations. Numerical calculations |35 , 36fl give support for a picture where these 



deviations are small up to the fragmentation regions, therefore we neglect them here. 

With the scaling ansatz for the longitudinal velocity we can relate the initial density 
and initial time to the final multiplicity distribution of hadrons by n^Tj = (dN/dr))/(7rR\). 
The temperature is therefore a function of 77, T« = Ti(i]), where we can identify the fluid 
and particle rapidities since the overall rapidity distribution is much broader than the 
thermal distribution at freeze-out. We parametrize the multiplicity distribution in the 
form 



dN ( dN\ . 9 lr% „, 

— = ( — ] exp(- V 2 /2a 2 ) , (11) 



drj \ drj 



where (dN/dr])o is the total multiplicity at rj = and a measures the width of the hadron 
rapidity distribution. As before, [jXJ] we assume the following initial conditions: at RHIC 
energies, \rj\ < 6 and we choose (dN/drj) Q = 2000 while for the LHC \r]\ < 8 and we take 
(dN/di]) = 5000. The width is estimated from the total energy of the final particles. We 
find a ~ 3 at RHIC and at LHC, o ~ 5. These rapidity densities, while somewhat high, 
are consistent with the (optimistic) high initial temperatures taken here, as well as with 
Parton Cascade Model calculations |3~l| . Like the scaling properties of the longitudinal 
velocity, the Gaussian form of the multiplicity distribution must break down at rapidities 
close to the phase space limit. There, however, the densities become too small for thermal 
production to be important. 
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3.1 Thermal Dileptons 

In the quark-gluon plasma, lepton pairs are expected to be produced by the annihilation 
process qq — > 7* — > l + l~, similar to Drell-Yan production in the initial nucleon-nucleon 
interactions. However, in the plasma, the nuclear quark distributions are replaced by 
thermal distributions. Neglecting the transverse flow the dilepton emission rate is then 



37, 38 



d^dydM 2 = ^e W [-(rr^cosh(y-rj)/T)} , (12) 

where T is the temperature and r\ is the rapidity of the medium. In a plasma, Fq = Y^ e q, 
while in an equilibrium hadron gas, assuming tttt — > p — > is the dominant channel, 
Fb = ^tti a / {{w? p — M 2 ) 2 + m 2 T 2 ). The mixed phase has a fractional contribution from 
each with F M = f (r - 1) [f F Q + ((r - 2)/ + 2)F H }. 

In addition to ttti — > p — > other hadronic processes contribute to the dilepton rate. 
In the region above the p mass they dominate over the p contribution which is cut 



off by the form factor. The assumption that vector mesons saturate the e + e~ — > hadrons 
cross section indicates that at the phase transition temperature the dilepton rate in the 
hadron gas could be very similar to that from the quark-gluon plasma [fiOfl . Neglecting 



these contributions underestimates the hadronic emission in the mass range between the 
p and J/ip mesons. 

The rapidity distribution for a given dilepton mass is 



dN a 2 2 f 3F Q f y dyjrqfl 

-ttR a 6( Sl -s Q )— ^ —j- -r Pi{x)e U 



dMdy 2tt 3 a [ v w ' M 3 J-y cosher] - y) 
+ 6{ Sl - s H )M^F M f_ Y drj (1 + 1) 



(13) 



+ eis . _ Sd fjk r mt h t* h )* Pi(x)e - rH 



where P\(x) = x 4 + 5x 3 + 15x 2 + 30a; + 30, x = M cosh(r/ — y)/T, Sj is the initial entropy 
density and sq and sh are the entropy densities of the plasma and the hadron gas at 
the transition temperature T c . The invariant mass distribution is obtained by integrating 
eq. (13) over the rapidity interval — Y < y < Y. Likewise the transverse momentum 
distribution for mass M is 



dN a 2 2 f f Y dr]dy(TiT 3 f 



-Z\Zi 
I 2c 



t 2 r Y 

+ 0(3, - s H )^-F M J drjdye- mTCOsHri - y)/T (14) 



+ 9{ 8i - s dcc )3F H J_ 



Y dr]dy(T H T% 



3 \2 



-Z\Z H 
' z dcc 



-y m\, cosh [r\ — y) 

where P 2 (z) = z 5 + 5z 4 + 20z 3 + QOz 2 + 120z + 120 and z = m T cosh(?7 - y) /T. 
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As defaults, we assume a three-flavor plasma and a hadron gas of massless pions with 
T c = 200 MeV and T dcc = 140 MeV. Then 

90 * Si > SQ ° r Sl < SH 

f S Q + (1 - f )s H S H < Si< Sq 

where 7^ is the number of degrees of freedom with jh = 3 and 7q = 16 + 21n//2 and 
rii = Si/3.6. The beginning of the mixed phase occurs at 

Tm I 1 fin T i = T c . 

The initial plasma content of the mixed phase is 

1 T>T C 
fo = { (si - s H )/ {s Q - s H ) Ti = T c . 
Ti < T c 

The beginning of the hadron phase is then 



th 



t r T > T 

'in 1 % _ c 

1 f m T < T r 



where r = sq/ sh = Jq/jh- There is no thermal contribution if Ti < Td ec - 

In our previous work, we used three different hypotheses to fix T and Ti [0. Since 
here we wish to determine if any plasma contribution is observable at all, we focus on the 
highest temperature scenario, Tj(3Tj) ~ he, yielding T ijmax ~ 515 at RHIC and 810 MeV at 
the LHC, similar to the results of Ref. |^TJ. The corresponding initial times are Tj ~ 0.2 fm 
and 0.08 fm respectively. Thermalization at such early times is a bold assumption, if taken 
literally. Following McLerran, Kapusta and Srivastava |5l| , we employ it to schematically 



characterize pre-equilibrium pair production |52|, 4"5fl . For our choice of initial conditions, 
plasma production dominates the thermal distributions over all observable phase space. 
However, other initial conditions with lower initial temperatures result in a much smaller 
thermal yield for M > 2 GeV. For example, if we choose instead T c = 150 MeV and 
^dec = 100 MeV, as recent lattice calculations with quarks indicate fH|, then the hadron 



gas contribution is considerably reduced. Since we have assumed a high initial temperature 
and study pairs with M > 2 GeV only, the hadron contribution is a small fraction of 
the yield. At M = 2 GeV, the yield is reduced by ~30% with T c = 150 MeV while 
higher masses are virtually unaffected. Reducing (dN/dr]) would reduce Tj, increasing 
the relative importance of the hadronic contribution and decreasing the plasma yield. 
Additionally, recent analyses of the initial conditions [[45], Q suggest that the quark 
density is too low for the quarks to be in chemical equilibrium. Therefore, even if the 
initial temperature is high, the dilepton yield would be significantly reduced because of 
the low quark density. 

We neglect the transverse expansion and describe the matter in the longitudinal direc- 
tion only. Transverse expansion mainly affects the later part of the mixed phase and the 
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hadron gas phase [47, 48|. If it is included, the contribution from the hadron gas phase 
is negligible for M > 2 GeV except at large rapidities where the initial density is too 
low to produce a plasma. Therefore a full three-dimensional hydrodynamic calculation 
would give a narrower thermal pair rapidity distribution than our longitudinal scaling 
approximation. 

We have so far discussed the properties of the lepton pairs. However, once acceptance 
questions are addressed, it is necessary to track each lepton separately. The space-time 
integration remains the same, but now we write more generally 



dN r d 3 k 1 d 3 k 2 d 3 p 1 d 3 p 2 .36\M\ 



2 



^ d A xd 3 pdy g dM 2 = J h)f{h) 16(2^ {h + h ~ Pl ~ V ^ 



6(M 2 - (p! + p 2 y)5(y g - y P1 + y P2 )S(p- n - n) , (15) 

where k\ and k 2 are the incoming quark and antiquark, p\ and p 2 are the outgoing leptons, 
and f(k) is the thermal (Boltzmann) distribution of the initial partons. Note that if the 
rapidity gap is integrated over, eq. (15) reduces to eq. (12). The px, mass, rapidity gap, 
and pair rapidity distributions for M > 2 GeV are shown in Figs. 8 and 9 for RHIC 
and LHC respectively. The number of thermal lepton pairs with M = 2, 4, and 6 GeV 
are given in Tables 1 and 2 for central collisions at RHIC and LHC. The rapidity gap 
distribution, dN/dy g dM does not increase with y g at fixed M as it does for the initial 
QQ production because the thermal partons are not aligned along a particular direction. 
However, the average px and rapidity of the thermal distributions follows the trend of the 
heavy quark decays: (py) increases with mass while (y) decreases. The average thermal 
pair mass and px are shown in Table 3 for M > 2 GeV at RHIC and LHC. 

3.2 Thermal Charm 

With a large initial temperature in the most optimistic scenario, Tj jmax ~ (1/3 — l/2)m c , 
significant thermal charm production may be expected, as previously suggested ||9|, pTi| ]. 
The thermal cc production rate can be found by making the replacement 

36\M q - q ^ u \ 2 -> l q \M q - q ^ 2 + j g \M gg ^ c - c \ 2 , (16) 

in eq. (15) where now p\ and p 2 are the four-momenta of the charm quarks and the 
quark and gluon degeneracy factors are 7 g = 3 x (2 x 3) 2 for three quark flavors and 
7 9 = (2 x 8) 2 /2 (the factor of 1/2 is needed to prevent double counting) [f)0|| . The 



matrix elements are given in [[5T|. We use m c = 1.2 GeV to be consistent with the initial 
cc production calculations and include only plasma productionP]. The fragmentation of 
the charmed quarks into DD pairs has also been incorporated, as in eq. (5). We find 
approximately one thermal DD pair at RHIC and 23 pairs in Pb+Pb collisions at LHC. 
These results are essentially in agreement with those of the ideal thermal case discussed in 



52] and [53]. The rate is high enough at the LHC for pa 500 lepton pairs to be produced 



6 In our previous paper 0, a normalization constant was left out of the calculation. The correction 
has the effect of reducing the thermal charm yield. Note also that a typographical factor of two is also 
missing from eq. (13) of that paper. 
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from uncorrelated decays. Again, an ideal background subtraction should remove the 
uncorrelated pairs. 

The thermal charm and resulting DD distributions are shown in Figs. 10 and 11 along 
with the dilepton yields. The number of thermal lepton pairs with M — 2, 4, and 6 
GeV are given in Tables 1 and 2 for central collisions at RHIC and LHC. As before, the 
number of pairs at each mass is averaged over mass bins of width 400 MeV. The average 
DD pair mass and px are somewhat larger than the corresponding cc averages. The DD 
pair mass is 1.5 GeV larger while (pt) is 400 MeV larger at RHIC and 600 MeV larger 
at LHC. These increases deplete high rapidity DD pairs with respect to cc production. 
The rapidity gap is also decreased, although not significantly. The average thermal pair 
mass and pt are given in Table 3 for both RHIC and LHC. The mass and distributions 
of lepton pairs resulting from correlated DD decays are similar in shape to the thermal 
dileptons and with a similar yield at the same lepton pair mass. However, the lepton 
pair rapidity distributions from DD decays are narrower. As in the case of initial cc 
production, the uncorrelated lepton pairs have a much broader mass distribution than 
the correlated pairs. 

4 Results and Conclusions 
4.1 Results 

We now compare the perturbative lepton pair production with thermal production. The 
mass distributions have already been shown in Fig. 1 for RHIC and LHC although without 
the uncorrelated DD decays. The pr and rapidity distributions for lepton pairs with 
M = 2, 4, and 6 GeV are shown in Fig. 12 for RHIC and 13 for LHC. Here the yield from 
the uncorrelated DD decays are included. The average lepton pair masses for M > 2 GeV 
and transverse momenta at both energies are given in Table 3. The average pr of the 
lepton pairs are given for M > 2 GeV only for Drell-Yan and thermal dilepton production. 
The average px from heavy quark decays, initial and thermal, are given for all M. 

Uncorrelated initial DD decays dominate the distributions, by several orders of mag- 
nitude at the LHC. However, the like-sign subtraction that removes it + it~ and KK decays 
from the continuum should also remove the uncorrelated DD decays. The correlated DD 
decays are part of the signal and will not be subtracted. These decays dominate the con- 
tinuum up to M = 10 GeV. At higher masses, the BB decays begin to be as important as 
correlated DD decays. The same results are observable in the pt and rapidity distribu- 
tions. At RHIC, the contributions from the initial hard processes are above the thermal 
contributions over all phase space except at M = 2 GeV and y > 4.5, as seen in Fig. 
12(b). The thermal dilepton and thermal charm contributions are somewhat above those 
of the BB and Drell-Yan pairs for M < 3 GeV at RHIC. Both thermal contributions have 
very similar distributions and yields at RHIC energies. The thermal dilepton and thermal 
DD distributions are also similar at the LHC although the thermal charm yield is larger 
than the thermal lepton pair yield. At the LHC, the BB decays produce more lepton 
pairs than the thermals for M > 2 GeV. However, the correlated thermal charm yield is 
above the BB decays at M = 2 GeV. The uncorrelated thermal charm yield is larger than 
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the BB decay rate. If all the heavy quark decays could be subtracted, there might be a 
small window of opportunity to observe thermal dileptons, both prompt and from charm 
decays, over Drell-Yan production at low pt, as seen in Fig. 12. Because RHIC is at a 
significantly lower energy than the LHC, the slopes of the rapidity and pt distributions 
of each of the contributions are somewhat different, even in the central region, as seen in 
Fig. 12. The difference in the slopes could perhaps help disentangle the dilepton sources, 
if very large rapidities could be measured. All the LHC pr and rapidity spectra have 
similar slopes, see Fig. 13, making differentiation more difficult. 

Our results suggest that it is very unlikely that the thermal qq annihilation signal 
can be extracted. However, the experimental acceptance has not been included. The 
acceptance should be smaller for lepton pairs from DD decays than for Drell-Yan or 
thermal dileptons. Since, especially for DD decays, large lepton pair mass implies a large 
rapidity gap, at least one of the DD decay leptons may be outside the finite detector 
acceptance. In particular, relatively few high mass lepton pairs from uncorrelated DD 
decays will be detected, significantly reducing the uncorrelated yield even before like-sign 
subtraction. The BB decay pairs will have a larger acceptance due to the increase in pair 
Pt over the DD — > decays at the same mass which reduces the rapidity gap. Therefore 
a judicious choice of kinematic cuts can reduce the initial DD acceptance relative to other 
dilepton sources that produce nearly equal but opposite transverse momentum leptons 
with pt < M. To show how the finite detector acceptance changes the yield, we choose 
some realistic cases to examine in more detail. 



4.2 RHIC: PHENIX 

The RHIC PHENIX detector is specially designed to measure electromagnetic probes. It 
consists of two central electron arms with a rather small acceptance, \t}\ < 0.35 and ±30° < 
cf) < ±120°, one muon arm in the forward direction [|53J] covering the pseudorapidity region 
1.1 < i] < 2.4, and a second muon arm in the backward direction with a similar angular 



coverage [55]. We assume that electrons with momentum greater than 1 GeV and muons 
with momentum greater than 2 GeV can be efficiently detected. Although we choose 
a rather high momentum cut-off, lower momentum leptons will be measured as well. 
Additional coverage from the electron arms together with the muon arm will partially fill 
in the rapidity gap between the two detector systems through electron-muon coincidence 
studies. The rapidity gap between the lepton pairs becomes particularly important when 
finite acceptance cuts are applied. Pairs with a small rapidity gap are more likely to be 
accepted, particularly in the central electron arms. 

The detector acceptance can have a substantial effect on the relative rates, as shown 
in Fig. 14 for (a) the central electron arms, (b) the forward muon arm, and (c) the 
combined e/i coverage for PHENIX. Note that some of the distributions are limited by 
statistics. For each system, the accepted lepton pair mass distributions are shown. The 
percentage of accepted pairs with M > 2 GeV and the average lepton pair mass from 
all our sources in this mass range are given in Table 4. In principle all finite masses are 
accepted, as suggested from the HH decay distributions in Figs. 4-7. However, since the 
Drell-Yan and thermal dileptons calculations are most reliable for M > 2 GeV, we use 
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M = 2 GeV as a lower bound on the accepted masses. Note that many of the DD decay- 
pairs, both correlated and uncorrected, have an invariant mass less than 2 GeV so that 
the additional mass cut to compare the number of accepted pairs on an equal footing 
significantly reduces the total acceptance, especially for the uncorrelated DD pairs. 

The finite acceptance strongly reduces the uncorrelated DD production relative to 
the correlated production, as seen in Fig. 14. In the central detector, DD decays have 
the smallest acceptance due to the relatively large rapidity gap between the leptons. 
Particularly, the high mass uncorrelated pairs are removed from the spectrum. The 
BB decays have the largest acceptance here because the combination of the relatively 
small pair rapidity and rapidity gap favors their detection. They will have the largest 
contribution to the continuum for M > 6 GeV, after the correlated DD signal is negligible. 

In the forward muon arm, the Drell-Yan and BB decays have very similar yields for 
M > 3 GeV since the BB acceptance will be decreased relative to the tightly correlated 
Drell-Yan and thermal production with their broader rapidity distributions. Although 
both correlated and uncorrelated DD decays are reduced relative to the smaller rapidity 
gap of the BB decays, the uncorrelated DD acceptance is reduced still further because 
of the larger rapidity gap between the uncorrelated lepton pairs. In general, the broader 
rapidity coverage in the muon arm increases the acceptance of both the correlated and 
uncorrelated DD pairs. After like-sign subtraction, the initial DD decays will dominate 
the spectrum for most of the pair masses studied. 

In Fig. 14(c) we show only the acceptance for correlated and uncorrelated initial DD 
decays to e/x pairs. The increased rapidity coverage of the combined system results in the 
acceptance of more low mass uncorrelated DD decays. We have not included the yield 
from any of the other sources since e/i pairs cannot be produced by correlated Drell-Yan 
and thermal production. The rate must be large enough for at least two lepton pairs to 
be produced per event for uncorrelated e/i pairs to be important (except for BB decays) 
and the yields from the other sources are small enough for such production to be unlikely 
(see Table 1). 

If the uncorrelated DD decays can be completely removed by a like-sign subtraction, 
then e/i coincidence is a good way to extract the correlated DD yield, which constitutes 
the charm signature. Even with a complete charm measurement, thermal sources will be 
hard to detect, especially the thermal dileptons. In the muon arm, Drell-Yan and BB 
decays have nearly the same rate and could be hard to separate. However, in the electron 
arms, for 6 < M < 10 GeV, BB decays are most important and could be removed from 
the lower mass continuum by a comparison of the slopes. The thermal DD signature may 
be measured below 3 GeV if the initial charm production can be reliably subtracted. At 
2 GeV in the central detectors, the accepted thermal DD yield is five times larger than 
the BB yield. At low mass and px, thermal charm could also be seen in the pair p T 
distributions. Measurements from all the detector systems must complement each other 
for reliable results. 
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4.3 LHC: ALICE 



At the LHC, only one heavy-ion detector, ALICE, is planned. It includes a central 
detector that will measure dielectrons with M < Mj/^ |56j and covering \rj\ < 0.9. We 
select leptons with momentum greater than 1 GeV. The ALICE collaboration has also 



proposed a forward muon spectrometer, with 2.4 < rj < 4 to cover higher mass pairs |57 
We consider only muons with momentum larger than 4 GeV in the muon arm. 

The mass distributions for (a) the central detector and (b) the forward muon spec- 
trometer are shown in Fig. 15. The relative rates are similar to the corresponding RHIC 
detector systems although, overall, the acceptances are larger in the ALICE detector, as 
shown in Table 5. In the central detector, this is probably due to the full azimuthal cover- 
age and the larger rapidity coverage. While the acceptance cuts substantially reduce the 
rate from uncorrelated DD decays compared to the correlated decays, the uncorrelated 
yield is still nearly an order of magnitude larger than the correlated yield for masses below 
4 GeV. This higher acceptance for uncorrelated pairs means than an accurate like-sign 
subtraction is crucial. An additional background comes from uncorrelated thermal DD 
decays. In the central detector, the uncorrelated thermal DD rate is as large as the 
initial correlated DD rate at low masses. However, in the muon arm, the BB decays 
are clearly the most important source of lepton pairs after the initial charm production 
due in part to the smaller rapidity gap between the leptons — it is more likely that both 
decay leptons will lie within the rapidity window of the detector. The acceptance does 
not decrease in the muon arm because the pair rapidity distribution is not significantly 
reduced with respect to the Drell-Yan and thermal pair distributions, as is the case at 
RHIC. The correlated thermal DD yield may also be observable at M ~ 2 GeV although 
the BB decay rate is within a factor of two here and more difficult to remove reliably. 
The average mass of the accepted pairs in ALICE are also given in Table 5. The observed 
trends are similar to the PHENIX results except that the conclusion that the thermal 
dilepton signal is unlikely to be measured is even stronger here although thermal charm 
may still be observable. 



4.4 Conclusions 

Charm production is the dominant source of dileptons in heavy-ion collisions, even with 
acceptance cuts, for M < 6 — 8 GeV. Uncertainties in QCD calculations may change the 
rates by a factor of two at RHIC and 3-4 at LHC, not enough to affect this conclusion. 
Charm is both signal and background because the multiple cc pair production results in 
substantial uncorrelated DD contributions to the background. We have only included 
heavy quark production by first collisions. However, multiple hard scattering in AB 
collisions can increase the charm yield before equilibration. 

Nuclear shadowing is not yet well understood for the gluon. If the shadowing effects 
can be mapped out in phase space, the detection of thermal signals could be improved. 
Since pA studies are planned at RHIC, nuclear shadowing could be mapped out in phase 
space. Such measurements are themselves important results. The effects are strongest 
at the LHC where the saturation of the shadowing curve is reached. Unfortunately this 
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saturation region is unlikely to be probed at RHIC and no corresponding pA measurements 
can be performed at the LHC so that the shadowing effects may be more difficult to 
interpret. It is clear that systematic studies of charm production in pp, pA and AB 
interactions at the same energy are needed to fully understand charm production. 

We stress that our work differers from previous efforts j|l|, 42] primarily in our esti- 
mate of the perturbative background. Kapusta, McLerran and Srivastava |41] assumed 
hydrodynamic initial conditions similar to ours, but concluded that thermal dileptons 
dominate the continuum below the T. We attribute this striking difference to their esti- 
mate of the initial hard scattering processes. First, they omitted the contribution from 
semileptonic charm decays. They also underestimated the Drell-Yan contribution by us- 
ing Duke-Owens parton distributions ||11|| , long obsolete. This same Drell-Yan estimate 
was also used in the comparison with dilepton production by the parton cascade model 

Our calculations show that charm cannot be omitted. The production of charm-decay 
and Drell-Yan dileptons in primary collisions in any parton cascade or hydrodynamic 
model must agree with perturbative QCD and, therefore, the rates should be similar 
to ours. Other initial charm calculations [[5^] are within a factor of two to three of 



ours, within the theoretical uncertainties involving the quark mass and scale. The cc 
production rate in the parton cascade model agrees with perturbative QCD in pp collisions 
but overestimates the importance of flavor excitation [TJ], leading to a larger charm yield 
than ours, particularly at RHIC However, we have omitted cascading of the c and 
c quarks in the high density medium, which can lead to energy loss |58| (similar to 'jet 
quenching' f59fl ). If this loss is sufficient, these quarks can be equilibrated with the flowing 
plasma. Since it is highly unlikely that all of the cc pairs can annihilate, cascading will 
not change the number of pairs appreciably. 

We expect that thermal charm will prove to be an experimentally accessible tem- 
perature probe at RHIC and LHC. On the other hand, we emphasize that thermal qq 
annihilation, perhaps the more familiar thermal signal, will be much more difficult to pick 
out. Thermal annihilation would be a more direct thermometer because the kinematics of 
the lepton pair specifies the off-shell photon's four momentum. However, the heavy quark 
and Drell-Yan contributions are too high for the steeply-falling thermal contributions to 
be extracted, unless the charm contributions can be reliably subtracted. 

What is the best way to measure charm? Coincidence measurements of efi can prove 
useful. Charm was first measured by this method at the ISR |60| and such coincidence 
measurements are planned for PHENIX Pairs of like-sign electrons may also offer 
a measure of uncorrelated charm production. Charm was measured with single electrons 
by a study of the e/ir ratio at the ISR [|61| . Additionally, semileptonic decays can be 
experimentally tagged and separated from direct production of lepton pairs with a vertex 
detector. If a detached vertex is observed for at least one of the leptons, then direct 
production of the pair can be ruled out. Such a vertex detector is planned for the STAR 
detector at RHIC. Another technique for reducing the signal levels from semileptonic 
decays is the use of selective kinematic cuts. Since the leptons from decays have a weaker 
correlation in rapidity or angle than those from directly produced pairs, cuts can be placed 
on these variables. While some signal events will be lost, the signal to background ratio 
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can be improved for large acceptance detectors. 
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A Appendix: Shadowing Effects on Initial Processes 

When the charged parton distributions are probed in deep-inelastic scattering with a 
nuclear target and compared to a deuterium target, the ratio Rp 2 = F^jF® has a char- 
acteristic shape as a function of x. The low x region, below x ~ 0.1, is referred to as 
the shadowing region, and between x ~ 0.3 and 0.7 is the EMC region. In both regions 
a depletion is observed in the heavy nucleus relative to deuterium. At very low x, Rf 2 



appears to saturate ||63| . Between the shadowing and EMC regions, an enhancement 
occurs, called antishadowing, where Rf 2 > 1- There is also an enhancement as x — > 1, 
assumed to be due to Fermi motion of the nucleons. The entire nuclear dependence is 
often referred to as shadowing. Although the behavior of Rp 2 is not well understood, 
the effect has been described by either an interplay of coherent and incoherent multiple 
scatterings in the target or a modification of the parton densities in nuclear matter. In 
any case, the effect can be modeled by an A dependent fit to the nuclear deep-inelastic 
scattering data and implemented by a modification of the parton distributions in the pro- 
ton. In this appendix, we show the effect of two different parameterizations of the nuclear 
parton densities to illustrate how the Drell-Yan and QQ distributions calculated for pp 
interactions might change in nuclear collisions at RHIC and LHC. 

In the central region at RHIC and LHC, the values of x probed are small enough for 
the hard processes we consider to be predominantly in the shadowing region. However, 
the momentum fractions increase with pair mass, transverse momentum, and rapidity. At 
y = and pt = 0, x ~ Mj y/s so that in the mass range 2 < M < 6 GeV, 0.01 < x < 0.03 
at RHIC and 3.6 x 10" 4 < x < 1.09 x 10~ 3 at the LHC. It follows that R F2 may change 
significantly at RHIC. At the LHC, x is small enough for the shadowing to be saturated 
over most measurable rapidities. The most important point to note is that the reduction 
in the Au+Au cross section is never significantly more than a factor of two for any of the 
processes considered. If the shadowing function can be mapped out in pA interactions at 
RHIC, as discussed in Ref. Q, then the corrections to AA interactions may be relatively 
clear for QQ and Drell-Yan production, especially since the lepton pairs from QQ decays 
should reflect the shape of the shadowing function in the same way as the QQ pairs 
themselves ||. 

The first parameterization is a general fit to the most recent nuclear deep-inelastic 
scattering data. The fit does not differentiate between quark, antiquark, and gluon mod- 
ifications and does not include evolution in Q 2 . It is not designed to satisfy the baryon 
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number or momentum sum rules. The functional form of Rp 2 is [15 



R 



F 2 



1 + 0.0134(1/2; - l/x sh ) 
a l + 0m27A^(l/x-l/x sh ) X<Xsh 

"emc b emc X -^sh X <L Xf erm ; j (1) 

/ 1 _ Xfcrmi \ 0.321 
R f\— — ) Xfermi < X < 1 



X 



where R s = a emc - b emc x sh , R f = a emc - frcmc^fermi, b emc = 0.525(1 - A~ 1/3 - 1.145,4~ 2/3 + 
0.93A- 1 + O.88A- 4 / 3 - 0.59A- 5 / 3 ), and a emc = 1 + b cmc x cmc . The fit fixes x sh = 0.15, 
iemc = 0.275 and £f erm i = 0.742. In Fig. 16(a) we show Rf 2 ( x ) f° r A = 197. In the 
nucleus, the parton densities are modified so that 

qf(x,fi) = R F2 (x)q p f (x, n) (2) 
g A (x,n) = R F2 (x)g p {x,n) , (3) 

where g/ represents both valence and sea quarks. Since Rf 2 ( x ) is scale independent and 
the parton densities are treated equivalently, the ratio of hard process production in 
Au+Au to pp collisions at the same energy is 

^ T ' M) = da pp /dyd PT dM " R *M R *M ■ W 

In Fig. 16(b)-(d) and Fig. 17 we show the nuclear effect on leading order calculations 
of heavy quark and Drell-Yan production. Figures 16(b)-(d) show R(pr) for single c and b 
mesons and R(M) and R{y) for DD and BB pairs at RHIC and LHC. At RHIC energies, 
the increase of R with px and M reflects the integration over the low xi, X2 midrapidity 
contributions as well as the growth of Rf 2 as x approaches the antishadowing region, 
reached at M ~ 20 GeV for DD and BB pairs. Note that R(M) is almost identical for 
DD and BB production, as it should be since at fixed M the same x values are probed. 
The B p T ratios are generally flatter because the change in R F2 (x) with p T is slower than 
the change in Rp 2 (x) for D production at the same energy. The ratio R{y) is nearly 
constant for DD production at RHIC, caused by the coincidence of x\ increasing toward 
the antishadowing region while X2 decreases into the saturation region. In contrast, the 
BB ratio decreases with rapidity at RHIC since at y = 2, x\ lies in the EMC region while 
x 2 ~ 0.001. Figure 17 shows R(y) and R{M) for Drell-Yan pairs at RHIC and LHC. At 
RHIC R(y) increases with mass at y = due to the increase of Rp 2 as x\ = x 2 grows 
from 0.01 to 0.03. At higher rapidities, Rf 2 {x2) moves toward the saturation region while 
Rp 2 (xi) passes through the EMC region and then rises again as x\ 1. Near the edge of 
phase space, the transition from the EMC region to the Fermi motion region can be seen. 
The phase space is depleted at higher rapidities for lower masses. At the LHC x\ remains 
in the shadowing region for y < 4. The ratios of rapidity integrated mass distributions 
rise slowly with mass, reflecting the broad rapidity distributions. 

The second parameterization modifies the valence and sea quark and gluon distribu- 
tions separately and includes evolution with the square of the momentum transfer, Q 2 , 



16] but is based on an older fit to the data using the Duke-Owens parton densities pl|. 
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The exact form of Rp 2 (x, Qq) is given in Ref. [16]. The initial scale is chosen to be Qq = 2 



GeV and the Q 2 evolution is done with both the standard Altarelli-Parisi evolution and 
with gluon recombination at high density. The gluon recombination terms do not strongly 
alter the evolution. In this case, the nuclear parton densities are modified so that 

q A (x,/i) = R v (x,/J l )q^(x,ii) (5) 
?sO,M) = R S (x, /i)q s (x, /J,) (6) 
g A (x,/j,) = R G (x,/j,)g p (x,iJ,) , (7) 

where qy = u v + d v is the valence quark density and qs = 2(u + d + s) is the total sea 
quark density and we assume that Ry and Rs affect the individual valence and sea quarks 
identically. The ratios were constrained in the model ]16] by assuming that Rp 2 ~ Ry 
at large x and Rp 2 ~ Rs at small x since xqy(x,fi) — > as x — > 0. We use ansatz 
1 for the gluons, Rp 2 ~ Rg f° r all x [IB|, since one might expect more shadowing for 



the sea quarks, generated from gluons, at small x. The parton densities satisfy baryon 
number conservation Jq 1 dx qy A (x, fx) = 3 and momentum conservation Jq 1 dx x(qy A (x, fi) + 
qs' A (%,fJ>) + g p ' A (%, /-0) = 1 at all /i. We have used the MRS D— ' densities with this 
parameterization instead of the original parton densities, leading to some small deviations 
in the momentum sum but the general trend is unchanged. 

Figure 18(a) shows the ratios Ry (solid curves), Rs (dashed curves), and Rq (dot- 
dashed curves) for the minimum, Q , and maximum, 10 GeV, values of Q for A = 200. 
Outside this range the ratios are fixed to those at 10 GeV. The valence quarks show little 
Q 2 evolution, the gluons the greatest. The sea quarks evolve more slowly than the gluons 
and, at Q = 10 GeV, Rs ~ Rf 2 °f the first parameterization, eq. (1). The ratios of AA 
to pp production of QQ and Drell-Yan pairs are shown in Fig. 18(b)-(d) and Fig. 19. The 
ratio R for this parameterization is not as straightforward to write down as in eq. (4) 
except for gg fusion. For example for qq — > QQ, 

R ^ Rs(xi)Rs(x 2 )^ q =u,d,s2qs(xi)q s (x 2 ) + [Ry(x 1 )Rs(x2)T, q=u4 q v (x 1 )q s (x 2 ) + (1 -> 2)] ^ 
^ q =u,d,s^Qs(xi)q s (x 2 ) + [^ q =u,dqv{xi)q s (x 2 ) + (1 -> 2)] 

where q s is a generic sea quark distribution and the \i dependence has been suppressed. 
The fast evolution of the gluons has the strongest effect on the charm and bottom pro- 
duction since QQ production by gluons is dominant. In fact, for bb production at RHIC 
energies, the shadowing effect has nearly vanished in R(pr) ■ The ratios for Au+Au to pp 
are shown in Fig. 18(b)-(d) as a function of quark p? and pair mass and rapidity. The 
trends are the same as for the first parameterization but, overall, the QQ distributions are 
not as strongly modified since the effect decreases for increasing p? and M. As seen in a 
comparison of Fig. 19 with Fig. 17, the shadowing effect is actually stronger for low mass 
Drell-Yan production with this parameterization due to the strong sea quark shadowing 
at Qq. 

In either case, the reduction in AA yield relative to pp due to shadowing is generally 
not larger than a factor of two, depending on the shadowing model. It is important to 
note that the total depletion is a dependent on both x and Q and is not a constant factor 
as a function of pt, M and y in either model. 
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Source 


M a = 2 GeV 


Mu = 4 GeV 


Mu = 6 GeV 


DD uncorr 


1.94xl0 _1 


4.5xl0" 2 


1.67xl0 -2 


DDcorr 


2.31xl(r 2 


2.10xl0- 3 


3.10xl0~ 4 


BB 


1.38xl0 -4 


6.85xl(r 5 


2.52xl(r 5 


DY 


5.29xl0~ 4 


3.86xl0~ 5 


8.06xl(r 6 


^th 


9.07xl(T 4 


7.6X10" 6 


1.4xl0" 7 


DD th 


9.29xl(r 4 


5.18xl(r 6 





Table 1: Number of lepton pairs per event from each of our sources in central Au+Au 
collisions at RHIC. 



25 







Source 


M u = 2 GeV 


M u = 4 GeV 


M u = 6 GeV 




4.80xl0 2 


1.90xl0 2 


1.07xl0 2 




1.40x10° 


1.69X10- 1 


2.95xl0~ 2 


BB 


2.05xl0~ 2 


1.15xl0~ 2 


4.98xl0" 3 


DY 


6.90xl0~ 3 


7.83xl0~ 4 


2.06xl0~ 4 




1.43xl0~ 2 


4.68xl0~ 4 


3.26xl0~ 5 


D -Dth ,corr 


4.8xl0" 2 


1.02xl0~ 3 


3.57xl0~ 5 


DDth ,uncorr 


1.42x10° 


4.32X10" 1 


2.03X10- 1 



Table 2: Number of lepton pairs per event from each of our sources in central Pb+Pb 
collisions at LHC. 
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200 GeV 


5.5 TeV 


Source 


(M u ) (GeV) 


M (GeV) 


(M u ) (GeV) 


(Pt,u) (GeV) 




3.73 


0.84 


4.57 


1.00 




2.41 


0.77 


3.03 


0.84 


BB 


3.95 


1.91 


4.20 


2.05 


DY 


2.45 


1.68 


2.68 


2.45 




2.44 


1.06 


2.61 


1.26 


-D-Dth.corr 


2.39 


0.85 


2.52 


0.98 


D -Dth,uncorr 






2.28 


0.97 



Table 3: Average lepton pair mass (for Mu > 2 GeV) and transverse momentum at RHIC 
and LHC. For the Drell-Yan and thermal lepton pairs, the average pr's are calculated for 
Mu > 2 GeV only while all the average pr's from heavy quark decays are for all masses. 
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e + e 






Source 


% Acc. 


(M ee ) (GeV) 


% Acc. 


(M m ) (GeV) 


% Acc. 


(M eM ) (GeV) 


DD uncorr 


0.01 


2.75 


0.07 


2.53 


0.084 


3.2 




0.032 


3.05 


0.42 


2.84 


0.08 


3.74 


BB 


0.60 


4.53 


1.7 


3.56 






DY 


0.16 


3.15 


3.7 


2.64 






^th 


0.23 


2.51 


3.3 


2.18 






DD th 


0.014 


2.47 


0.12 


2.58 







Table 4: Percentage of lepton pairs with Ma > 2 GeV accepted in the PHENIX detector 
and their average pair mass. We have included the central electron arms, the forward 
muon arm, and a combination of electrons and muons from both detectors. 
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e + e 






% Acc. 


(M ee ) (GeV) 


% Acc. 


(M m ) (GeV) 




0.093 


2.69 


0.12 


2.58 




0.44 


3.10 


0.67 


2.92 


BB 


4.42 


4.33 


3.58 


4.01 


DY 


3.6 


2.62 


5.27 


2.58 




4.9 


2.55 


4.76 


2.22 




0.46 


2.52 


0.29 


2.34 


DD^\y jU ncorr 


0.33 


2.70 


0.056 


2.49 



Table 5: Percentage of lepton pairs with M tt > 2 GeV accepted in the ALICE detector 
and their average pair mass. We have included the central detector and the proposed 
forward muon arm. 
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Figure Captions 



Figure 1. The lepton pair mass distributions are given in (a) for central Au+Au collisions 
at RHIC and (b) for central Pb+Pb collisions at LHC. The contributions are: Drell- 
Yan (dashed) and thermal dilepton (dot-dashed-dashed) production and thermal DD 
decays (dotted), as well as initial correlated (dot-dashed) DD and BB (dot-dot-dashed) 
production and decay. Note that lepton pairs from uncorrelated initial DD decays (as 
well as uncorrelated thermal DD decays at the LHC) have not been included here but 
are a very large contribution to the continuum. 

Figure 2. The NLO Drell-Yan pair production rate in central Au+Au collisions at y/s = 
200 GeV. The pair p T (a), mass (b), and rapidity (c) distributions are shown. The p T 
distributions are given at y = (solid) and y = 2 (dashed) for M — 2, 4, and 6 GeV (the 
upper, middle, and lower sets of curves respectively). In (c) the rapidity distributions are 
calculated for M = 2 GeV (solid), 4 GeV (dashed), and 6 GeV (dot-dashed). 

Figure 3. The same as Fig. 2 for central Pb+Pb collisions at y/s = 5.5 TeV. 

Figure 4. The rate of initial cc (solid) and DD (dashed) pair production to NLO in central 
Au+Au collisions at y/s = 200 GeV. The pair px (a), mass (b), rapidity gap (c), and pair 
rapidity (d) distributions are shown. Additionally, the lepton pair pt, mass, and rapidity 
distributions from correlated (dot-dashed) and uncorrelated (dotted) DD pair decays are 
shown for central Au+Au collisions. 

Figure 5. The same as Fig. 4 for central Pb+Pb collisions at y/s = 5.5 TeV. 

Figure 6. The same as Fig. 4 for the bb and BB production rate to NLO in central Au+Au 
collisions at y/s = 200 GeV. 

Figure 7. The same as Fig. 6 for central Pb+Pb collisions at y/s = 5.5 TeV. 

Figure 8. The thermal lepton pair rate in central Au+Au collisions for M > 2 GeV 
and y/s = 200 GeV. The lepton pair pt (a), mass (b), rapidity gap (c), and rapidity (d) 
distributions are shown. The px, rapidity gap, and rapidity distributions are given for 
M = 2 (solid), 4 (dashed), and 6 (dot-dashed) GeV. 

Figure 9. The same as Fig. 8 for central Pb+Pb collisions at y/s = 5.5 TeV. 

Figure 10. The thermal cc (solid) and DD (dashed) pair production rate in central Au+Au 
collisions at y/s = 200 GeV. The pair p T (a), mass (b), rapidity gap (c), and rapidity 
(d) distributions are given. The lepton pairs from thermal DD decays are shown in the 
dot-dashed curves. 
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Figure 11. The same as Fig. 10 for central Pb+Pb collisions at ^/s = 5.5 TeV. Uncorre- 
lated thermal DD decays are also included (dotted curves). 

Figure 12. The contributions to the dilepton spectrum in central Au+Au collisions at 
y^s = 200 GeV for pairs with M — 2, 4, and 6 GeV. The p T distributions are given 
in (a), (c), and (e) while the pair rapidity distributions are shown in (b), (d), and (f). 
The distributions from Drell-Yan (dashed) and thermal dilepton (solid) production and 
thermal DD decays (dotted), the initial correlated (dot-dashed) and uncorrelated (dot- 
dashed-dashed) DD and BB (dot-dot-dashed) production and decay are included. 

Figure 13. The same as Fig. 12 for central Pb+Pb collisions at y/s = 5.5 TeV. In addition 
to the contributions shown in Fig. 12, uncorrelated thermal DD decays are shown in the 
dot-dot-dashed-dashed curves. 

Figure 14. The mass distributions of pairs accepted into the PHENIX detector in central 
Au+Au collisions at RHIC. The acceptance cuts are shown for pairs in (a) the central 
electron detector, (b) the forward muon arm, and (c) pairs formed when an electron is 
accepted into one of the PHENIX central arms and an opposite sign muon is accepted into 
the forward arm. The contributions are: Drell-Yan (dashed) and thermal dilepton (solid) 
production and thermal DD (dotted), initial correlated (dot-dashed) and uncorrelated 
(dot-dashed-dashed) DD, and initial BB (dot-dot-dashed) production and decay. In (c), 
only the mass distributions for correlated and uncorrelated DD — > e ± /i^X decays are 
shown. 

Figure 15. The same as Fig. 14 for pairs accepted into the ALICE detector in central 
Pb+Pb collisions at LHC. The accepted distributions are given for (a) the central detector 
and (b) the proposed forward muon arm. In addition to the contributions shown in Fig. 
14, uncorrelated thermal DD decays are shown in the dot-dot-dashed-dashed curves. 

Figure 16. (a) The shadowing function Rf 2 (x) fc> r A — 197. We also show the charm 
and bottom Au+Au to pp production ratios as a function of (b) heavy quark px, (c) QQ 
invariant mass and (d) QQ pair rapidity. Charm production is given in the solid curve at 
RHIC and the dashed curve at LHC. Bottom production is shown in the dot-dashed curve 
for RHIC and the dotted curve for LHC. Both Au+Au and pp production is calculated 
at v/s = 200 GeV for RHIC and ^/s = 5.5 TeV at the LHC. 

Figure 17. The ratio of Drell-Yan rapidity and mass distributions in Au+Au to pp col- 
lisions at yfs = 200 GeV at RHIC and 5.5 TeV at the LHC. The ratio R(y) is given for 
M=2 (solid), 4 (dashed) and 6 (dot-dashed) GeV for RHIC (a) and LHC (c). The mass 
R(M) is also given for RHIC (b) and LHC (d). 

Figure 18. (a) The shadowing functions Ry(x) (solid), Rs(x) (dashed) and R g (x) (dot- 
dashed) for Q = 2 (lower curves) and 10 GeV (upper curves) for A = 200. We also 
show the charm and bottom Au+Au to pp production ratios as in Fig. 16 but with the 
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shadowing functions given in (a). 

Figure 19. The same as in Fig. 17 with the shadowing functions as in Fig. 18(a). 
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